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Abstract 

We have measured quantum interference between two single microwave photons 
trapped in a superconducting resonator, whose frequencies are initially about 6 
GHz apart. We accomplish this by use of a parametric frequency conversion process 
that mixes the mode currents of two cavity harmonics through a superconducting 
quantum interference device, and demonstrate that a two-photon entanglement 
operation can be performed with high fidelity. 

Two identical optical photons simultaneously sent through the two input ports of a 
semi-transparent beam splitter will exit together through the same output port. This is 
a consequence of a complete destructive quantum interference between the two possible 
paths in which the photons leave by different outputs. This effect was first experimentally 
observed in parametric down-conversion [1, 2], and more recently by use of various single- 
photon sources [3, 4, 5]. It is at the heart of dual-rail linear optical quantum computing 
schemes [6, 7], where information is encoded on different spatial modes. However, these 
schemes can be extended to any type of system composed of linearly-coupled bosonic 
modes. 

Parametric frequency conversion (PFC) enables the creation of a tunable direct cou- 
pling between quantized linear modes of different frequencies via a time-varying, or a 
nonlinear element pumped at their frequency difference. It realizes an active beam split- 
ter (BS) component between these modes of different frequencies, with an experimentally 
tunable mixing angle [8, 9, 11, 10]. As such, it has been proposed as a way to imple- 
ment interferences [12, 13] between photonic quantum bits (qubits) encoded in frequency 
modes. In optics, PFC can be induced by intensely pumping a nonlinear crystal, and 
has been used to transduce a quantum state of light from one wavelength to another 
[14], down to the single-photon level [15, 16]. Superconducting circuits are attractive sys- 
tems to manipulate photons by means of parametric processes based on the microwave 
modulation of a Josephson element. The dynamics of the parametric coherent exchange 
of a single excitation between two superconducting qubits [17], or two linear resonators 
[10], far-detuned in frequency, has been recently measured. Parametric interactions can 
be combined with the tools developed for circuit quantum electrodynamics (cQED) [18] 
that enable the preparation and the read out of the quantum number states stored in 



a resonator using a superconducting qubit [19, 20]. Here, we prepare two single pho- 
tons of different frequencies in two resonant modes of a superconducting cavity that are 
parametrically coupled by a flux-modulated superconducting quantum interference de- 
vice (SQUID). We demonstrate quantum interference between these two single photons 
by implementing a Mach-Zehnder interferometer based on this parametric mode-mixing. 
Finally, we characterize the efficiency of the photon-pair preparation by measuring the 
fidelity of the two-photon entangled state created after interference on a semi-transparent 
(50:50) BS, and also give, as a comparison, the results obtained for single-photon inter- 
ference. 

Our circuit consists of a A/4 coplanar stepped- impedance resonator terminated to 
ground through a SQUID (Fig. la and lb) [10]. The other end is capacitively coupled to 
a 50 Q transmission line to allow spectroscopic measurements of the cavity by microwave 
reflectometry. At this node, the cavity is also coupled with a strength g q x = 17 MHz 
to a flux-biased phase qubit that is used to load (and detect) microwave photon number 
states in (from) one mode of the cavity [20]. An inductively coupled coil on-chip, referred 
as the "pump", allows the application through the SQUID loop of a global flux <3> sq = 
$gq + $^(t) composed of a small microwave modulation 3>^(t) on top of a d.c. bias &^q- 
The SQUID acts as a flux-dependent lumped-element inductor L sq ($gq) that slightly 
modifies one of the boundary conditions of the cavity [21, 22, 23]. The resulting flux 
dependence of the resonance frequencies i>o and z/i of the first two modes and 1 is 
shown in Fig. lb. The experiments reported here are performed at a flux bias $^ = 
0.21$ (point "A" in Fig. lb) where v£ = 3.8296 GHz and = 9.6087 GHz. The 
extracted loaded quality factors are about 16000 for both modes, high enough to perform 
dynamical manipulation of photon states trapped in the cavity via PFC [10]. By applying 
a small microwave drive $^q*(t) = 5$ cos(27r^ p t + ip p ) at a resonant pump frequency 
u p = Ui — Uq , with a pump phase ip p and an amplitude 5$, PFC is induced between 
the modes and 1, with a parametric coupling rate gio proportional to <5$. We apply 
a (J$ ~ 0.019 $o fhrx modulation, which yields to gio = 11.3 MHz. At a quantum 
level, PFC is described in terms of annihilation and creation operators ao,i and a\ 1 
of quantized microwave excitations, namely photons, of these harmonic modes, by the 
following interaction Hamiltonian [8] : 

Hj = hg 10 {aiaie iv » + a a\ e - iiPp ) (1) 

in the doubly rotating frame defined by the unitary transformation e 2j7r ^° a o a ° +Vl 
Hi describes a BS operation between modes and 1, with an effective mixing angle 
9 P = 47rgioAr p that can be modulated by varying the duration Ar p of the parametric 
drive. We consider the evolution of the | Ooli ) and |loli) states (the subscript denotes 
the mode number) under the unitary operation U(6 p ,(p p ) associated with the parametric 
interaction Hj, for 6 P G [0;27r]. The 6 p = n/2 operation is analogous to a 50:50 BS 
and creates the maximally entangled photon number state, a so-called NOON state [24], 
with one (N=l) or two (N=2) photons |^n) = [ ~ e iN{ ^ +n ^ |N 0i) + |0 o Ni)], 

where /n(<£ p ) is a global phase. In the N=2 case, destructive two-photon interference 
leads to the disappearance of | l li) in the final state, and the creation of \ip2) illustrates 
the bosonic nature of the photons, which ideally exit together through either "frequency 
output port". Note that the preparation of NOON states between two superconducting 
resonators has been demonstrated in ref. [25]. The method used in their work is different, 
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as it relies on transferring N times the single-excitation entanglement between two phase 
qubits into two distinct microwave resonators, the final state being characterized by a 
bipartite Wigner tomography. 

We prepare the initial state of the cavity by use of a phase qubit [26, 20, 19] (Fig. 
2a). The qubit frequency u q can be tuned from 7.5 GHz to 12.5 GHz by means of an 
external flux $ q . The qubit state is manipulated with microwave pulses and measured 
with an inductively coupled d.c. SQUID. By applying a fast flux pulse to the qubit, the 
excited state |e) state preferentially tunnels to a different flux state, which is indicated 
by a shift in the d.c. SQUID critical current. When the qubit is in the ground state, 
the tunneling probability P t is about 0.06; this probability increases linearly with the |e) 
state occupancy, up to 0.8 when the qubit is prepared in its excited state. We excite the 
qubit far off resonance from mode 1 (at z/° ff = 8.775 GHz), then bring it into resonance 
in order to swap the qubit excitation to the cavity, which prepares |Ooli). To prepare 
the | loli) state, we transfer the photon in mode 1 to mode by applying a 9 P = n 
pump pulse, then prepare and transfer another excitation from the qubit to mode 1. In 
practice, the relaxation of the qubit and of the cavity prevent us from preparing these 
initial pure Fock states. We prepare instead a mixed state described by the diagonal 
terms of an initial density matrix denoted p twm (9 p = 0) (p smgle (# p = 0)) when targeting 
the | loli) ( 1 0o 1 1 ) ) initial state, and expressed in the |m ni), m + ni < 2 (m + ni < 1), 
photon Fock states basis. We analyze the state of mode 1 after the U (9 P , (p p ) evolution 
using the qubit as probe. After applying the (9 p , <p p ) pump pulse, we bring the qubit 
into resonance with mode 1 for a varying interaction duration Ar ql and measure the 
qubit tunneling probability P t (Ar q i) (Fig. 2b and 2c). The probabilities P n (9 p ,(p p ) to 
measure n photons in mode 1 are given by the weights of the components scaling as 
-\/ng q i in the single-photon oscillations between the qubit and the mode 1 that interact 
during Ar ql (The P n (9 p , ip p ) are extracted by simultaneous least-squares fits to P t (Ar q i) 
with a Lindblad-type equation). They are related to the probabilities p mn (6 l p, <p p ) = 
(m iii \ p(9 p ,(p p ) |m n!) to simultaneously measure m and n photons in the mode and 
1 following P n (9 p ,ip p ) = J2 m +n<N Pmn(9 p ,(p p ), with N=l (N=2) in the single-photon 
(twin-photon) case. The extracted P n (# P , 0) (the phase of the first pump pulse is set to 
by convention and will not be specified in most cases) are shown in Fig. 2b and 2c. 
For both the single-photon and the twin-photon case, we measure P\{2-k) » 0.92 Pi(0) 
and can reasonably neglect the effect of decoherence over a pump period to simplify 
the analysis. We thus infer the photon number probability of the mode by use of a 
pump swap operation, and consider that Pno(#p) = Pon(# p + 7r). The measured initial 
population poi(0) and the estimated initial population pn(0) are shown in Table 1. In 
particular, when starting with p twm (0), we observe a \/2g q i component in the vacuum 
Rabi oscillations after a tt/2 (or a 3ir/2) pump pulse with a weight indicating that we 
have prepared at least Po2(tt/2) +£>2o( 7r /2) = 0.45 ± 0.04 of the | loli) state in the initial 
mixture. We characterize the quality of the NOON state created after a tt/2 pump 
pulse by performing Ramsey interferences [27, 28] (Fig. 2d). After preparing p smgle (|) 
or p twm (|), we apply a second n/2 pulse differing by a tunable phase A<p p .We measure 
with the qubit the final probability P n (|,0; |,A</? P ) to have n photons in mode 1 after 
this Ramsey sequence, as a function of A(p p (Fig. 2e and 2f). The coherent superposition 
of |OoNi) and |NoOi) induces a sinusoidal component oscillating at NA<p p in the Ramsey 
fringes pattern. More precisely, the consecutive application of C/(|,0j and [/(|,A<p p J 
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Table 1: Characterization of single and two-photon states 





Poi 


Pio 




Coil 


p sin § lc (0) 
p singlc (§) 


0.79(2) 
0.38(2) 


0.03(2) 
0.36(2) 




0.37(2) 




P02 


P20 


Pn 


C 02 


p twin (0) 
p twin (f) 


< 0.01 
0.22(2) 


< 0.01 
0.23(2) 


^ 0.45(4) a 


0.21(2) 



a This value is a lower boundary of Pn(0) given by the 
sum of the p2o( 7r /2) and po2(n/2). 
b The number in parentheses is the uncertainty on the 
last digit. 



projects the coherences con = (OoNi|p(|) |N 0i) = |coN|e lv3 ° N onto the final measured 
probability expressed by 

Pl (l' ' \ A ^ p ) = ~l c oi| sin ( A( ^p _ Voi) (2) 

+ |c 02 | cos(2Av9 p - ip 02 ) + Ai, 

where A\ is a constant depending on the p m n(§)> an d with c 02 = in the single-photon 
case. The main results are summarized in Table 1. From the two types of measurements 
shown in Fig. 2, we can estimate the fidelity J-"n = (V^lpd ) I^n) = ^(pon +Pno) + |con| 
of the NOON states preparation with our protocol combining transfer of an excitation 
from the qubit to the cavity, and mode coupling with PFC. We find T\ = 0.74 ± 0.04 
and Ti = 0.43 ± 0.04. The fidelity of the NOON state preparation is limited by the 
probability of initially preparing a single photon in mode 1 (poi(0)), or a single photon 
in each mode (pn(0)). The fidelity T\ is still higher than 0.5 by a significant amount, 
which is enough to claim deterministic entanglement between the two modes via a single 
photon [27]. Nevertheless, we can estimate the fidelity J 7 ^ of the N-photon entanglement 
operation, defined as the conditional probability J 7 ^ = ^n/ S m+n= NPmn(0), to create a 
NOON state, knowing that we started from one (N=l case) or two (N=2) photons in 
the cavity modes (Note that in the original Hong-Ou-Mandel experiment, this type of 
renormalization is achieved by working in the coincidence basis [1]). Due to relaxation of 
the qubit and the cavity during the pulse sequence, we have about 20 % probability to lose 
the excitation when transferring it from the qubit into mode 1 (Table 1). This leads to 
T\ = 0.90 ±0.06. In the N=2 case, we can place an upper bound on J2 m+n=2 Pmn{Q) using 
the probability Poi(0) = 0-62 ± 0.04 to load sequentially two excitations from the qubit 
to mode 1. This implies that T 2 is higher than 0.7 ± 0.06, which unambiguously proves 
the existence of two-photon quantum correlations in the created state. A more realistic 
upper bound of about 0.55 can be placed on X] m +n=2Pmn(0) by taking into account the 
loss during the transfer to the mode of the first photon loaded in mode 1, as well as its 
decay while preparing the second photon in mode 1. 

We also measure the decay of the NOON states by varying either the delay after a 
7r/2 pulse in order to access the decay of the populations Pon(|) and Pno(§ ) , or the delay 
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between the two 7r/2 pump pulses of the Ramsey interferometer in order to access the 
decay of the coherences con- We find that coi and C02 decay exponentially with respective 
characteristic times r| C[)1 | « 370 ns and T| C()2| « 190 ns, close to that of the harmonic mean 
2(1/t Pno + 1/t Pon ) _1 of the relaxation times t Pno and r Pm of the populations (see Fig. 
3), indicating that decoherence is limited by relaxation. We also observe, as expected 
[29], a decay rate of the two-photon entanglement approximately twice as fast as the 
single-photon one. 

To conclude, we have demonstrated quantum interference between two single pho- 
tons of significantly different microwave frequencies by use of PFC. The contrast of the 
observed interference is limited by the fidelity of the initial twin-photon state prepara- 
tion. Nevertheless, the fidelity of the two-photon entanglement operation is higher than 
70 %. Improving the relaxation time of the qubit and of the cavity would increase the 
quality of the two-photon entanglement calibration, as well as the photon-pair creation 
efficiency This could make this type of parametric device suitable for processing quantum 
information with photonic qubits encoded in microwave resonator modes. 
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Figure 1: (a) A A/4 CPW resonator with a boundary condition that is varied by a 
SQUID is coupled to a phase qubit. The flux through the SQUID is d.c.-biased and 
microwave modulated via an inductively coupled pump line. The cavity is coupled to 
a 50 Q transmission line through a small capacitance C c ~ 1.8 fF to allow reflectome- 
try measurements. We manipulate the qubit state with microwave pulses and tune its 
frequency with a flux bias line. Measurements are performed in a dilution refrigerator 
operated at 35 mK. (c) Measured phase Arg(Sii) of the reflection coefficient on C c as a 
function of both d.c.-flux in the SQUID and probe frequency. 
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Figure 2: (a) Microwave pulse sequence and optical analog for the states' evolution under 
PFC. Qubit, mode and mode 1 are respectively represented in green, red and blue. The 
mode mixing via the parametric pump is depicted in purple. (b, c) On the left, 2D plots 
of the tunneling probability P t (Ar q i,9 p ) when starting with p smgle (0) (b) or p twin (0) (c). 
On the right, the extracted probabilities P n to measure n photons in mode 1 are plotted 
as a function of 9 p . The error bars represent both the switching measurement statistical 
error and standard deviations from fits to the data, (d), Microwave pulse sequence and 
optical analog for the NOON states characterization. We apply a Ramsey sequence, i.e. 
two ix 1 1 pump pulses differing by a phase A<p p , on the initial state, (e, f), Ramsey fringes 
on the p single (0) (e) and on the p twin (0) (f) states. On the left, 2D plots of the tunneling 
probability P t (Ar q i, A<p p ). On the right, the extracted probabilities P n are plotted as 
a function of A<p p , and fitted (solid curves) by a cosine-sine decomposition restricted to 
frequencies of 0, 1, ...N. 
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Figure 3: Decay of the single and two-photon NOON states. The amplitudes of the 
coherence con, of the probabilities Pon(§) and PNo(f), are respectively plotted in dots, 
squares and diamonds, in logarithmic scale, as a function of the delay between either 
the two 7r/2 pulses of the Ramsey sequence for the coherence decay measurement, or the 
delay after a simple n/2 pulse for the populations decay measurements. These data are 
well fitted (solid curves) by an exponential decay. 
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